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ABSTRACT 

We study the temperature structure of the intergalactic medium (IGM) using a large cosmological N- 
body/SPH simulation. We employ a two-temperature model for the thermal evolution of the ionized gas, 
in which we include explicitly the relaxation process between electrons and ions. In the diffuse, hot IGM, 
the relaxation time is comparable to the age of the Universe and hence the electron temperature in post- 
shock regions remains significantly smaller than the ion temperature. We show that, at the present epoch, 
a large fraction of the warm/hot intergalactic medium (WHIM) has a well-developed two-temperature 
structure, with typical temperature differences of order a few. Consequently, the fraction of metals in 
various ionization states such as Ovi, Ovu, and Oviii, as well as their line emissivities, can differ locally 
by more than an order of magnitude from those computed with a single-temperature model, especially in 
gas with T ~ 10 7 K. It is thus necessary to follow the evolution of the electron temperature explicitly to 
determine absorption and emission by the WHIM. Although equipartition is nearly achieved in the denser 
intracluster medium (ICM), we find an appreciable systematic deviation between the gas- mass weighted 
electron temperature and the mean temperature even at half the virial radii of clusters. There is thus a 
reservoir of warm (T e < 1 keV) gas in and around massive clusters. Our results imply that relaxation 
processes need to be considered in describing and interpreting observational data from existing X-ray 
telescopes as well as from future missions designed to detect the WHIM, such as the Diffuse Intergalactic 
Oxygen Surveyor and the Missing Baryon Explorer. 

Subject headings: cosmology: theory - intergalactic medium 



1. INTRODUCTION 

The distribution of baryons in the Universe remains one 
of the puzzling issues in modern cosmology. At the present 
epoch, the total amount of baryons inferred from a cen- 
sus of Hi absorption, gas and stars in galaxies, and X-ray 
emission from hot gas in galaxy clusters is far smaller than 
that predicted by nucleosynthesis calculations (Persic & 
Salucci 1992; Fukugita, Hogan & Peebles 1998) and by 
measurements of the cosmic microwave background radi- 
ation (Spergel et al. 2003). Hence, it is believed that a 
large fraction of the baryons lies in an as yet undiscovered 
dark state. 

Three dimensional hydrodynamic simulations of cosmic 
structure formation (e.g., Cen & Ostriker 1999; Dave et 
al. 2001; Croft et al. 2002) predict that about 30-40% of 
the baryons reside in the so-called warm-hot intergalactic 
medium (WHIM) at the present epoch. This gas is mostly 
shock-heated to a temperature of ~ 10 5 -10 7 K during large 
scale structure formation, and this relatively low temper- 
ature makes its thermal emission difficult to detect with 
conventional X-ray probes. 

A variety of observational approaches have been consid- 
ered for studying the WHIM, including both ultraviolet 
(e.g., Furlanetto et al. 2003, 2004) and X-ray lines. With 
respect to the latter possibility, there have been several 
tentative claims that the WHIM has been detected locally 
in absorption (e.g., Nicastro et al. 2002) and in emis- 
sion (Finoguenov et al. 2003). Moreover, Yoshikawa et al. 
(2004) have proposed using Ovii/Ovin lines to examine 



the WHIM in detail with high spectral resolution X-ray 
detectors. They argue that, under a variety of assump- 
tions, about half of the WHIM (by mass) can be detected 
via oxygen line emission (see also Fang et al. 2004). While 
the proposed missions appear promising, the detectability 
depends crucially on the assumed metallicity and the rel- 
ative abundances of particular ions. 

Previous theoretical studies used a single-temperature 
model for the WHIM, under the assumption that equipar- 
tition is fully achieved. However, the outskirts of galaxy 
clusters may actually have a two-temperature structure 
(e.g., Fox & Loeb 1997; Takizawa 1998; Courty & Alimi 
2004). Near galaxy clusters, infalling gas shock- heats to 
roughly the virial temperature as kinetic energy is con- 
verted into thermal energy. Because the infalling ions 
carry nearly all of the bulk kinetic energy (exceeding that 
contained in electrons by a factor ~ m p /m e ), the post- 
shock electrons can gain thermal energy only via collisions 
with ions. If collisional relaxation takes a time compa- 
rable to or longer than the age of the Universe, the two 
temperatures, i.e. those of the electrons and ions, remain 
different. One can easily show that the relaxation time is 
rather long for diffuse warm-hot gas around clusters and in 
large-scale filamentary structures. Since the relative abun- 
dance of metal ions and their emissivities are sensitive to 
the electron temperature, this has important implications 
for future observations of the WHIM, particularly for the 
hotter component. 

In this Letter, we study the thermal evolution of the 
intergalactic medium (IGM) using a large hydrodynamic 
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Fig. 1. — The distribution of baryons (left), the warm/hot component with 10 6 < T < 10 7 K (middle), and gas with T c < 0.57] (right) in 
a slab of 100 X 100 X 20 (/i _1 Mpc) 3 at z = 0. The vertical bar in the middle panel indicates a 'super-cluster' region, for which we compute 
the soft X-ray intensity (see section 4 and Figure 5). 



simulation of structure formation. We show that a consid- 
erable fraction of the WHIM indeed has a two-temperature 
structure, particularly in and around rich clusters. In these 
regions, the electron temperature is typically smaller than 
the ion temperature by a factor of a few, substantially 
modifying estimates of the abundances of ions and their 
emissivities. 

2. THE Af-BODY/SPH SIMULATIONS 

For our simulations, we use the parallel Tree-PM/SPH 
code GADGET2 (V. Springel, in preparation), which em- 
ploys a fully conservative scheme for integrating the equa- 
tions of motion (Springel & Hernquist 2002). We imple- 
ment non-equipartition processes between ions and elec- 
trons following Fox & Loeb (1997) and Takizawa (1998). 
Equilibrium in an electron-proton plasma is achieved in 
the following manner. After passing through a shock, elec- 
trons and ions thermalize into (separate) Maxwellian dis- 
tributions on equilibration timescales t ce , ta, with ta/t ee ~ 
(mj/rTio) 1 / 2 ~ 43 (Spitzer 1962). Here, m c and m,i are 
the electron and ion masses and we have assumed pro- 
tons dominate the ionic component. Equipartition be- 
tween protons and electrons is achieved on an even longer 
timcscale t e i ~ (rrii/m e )t ee . Hence, we assume that elec- 
trons and ions quickly achieve Maxwellian distributions 
with temperatures T c and Ti, respectively, and consider 
only the non-equipartition effect. Spitzer (1962) showed 
that the appropriate timescale is given by 

3m c mi f kT c kTi \ 3 ^ 2 

tci= 8(2 7 r) 1 /2 n ^2 e 4 lnA \~^ + ^~ ) ' (1) 

where e is the electric charge, Z\ is the charge of an ion, 
7i; is the ion number density, and In A is the Coulomb log- 
arithm, which is given by 

lnA = 37.8 + ln(-5-^ --bx( T V (2) 

Vl0 7 K/ 2 V10- 5 cm- 3 / V ; 

For a hot, fully ionized plasma, 

t ei = 6.3 x 10V - 10 _ 7 3 K rA MOV < 3 > 

(n/10 a cm ■ 3 )(lnA/40) 



and this timescale can be comparable to the Hubble time 
in regions with T ;> 10 K and overdensities 10 — 100 with 
respect to the cosmic mean. The evolution of the electron 
temperature is given by 

ltT = ^— + {l - l) ^Tf (4) 

We do not consider heat conduction on the assumption 
that it will be suppressed by tangled magnetic fields. Some 
other electron heating processes driven by plasma insta- 
bilities have also been proposed (Laming 2000), but de- 
tailed calculations in the context of supernova remnants 
suggest that these mechanisms convert less than 20% of 
the bulk kinetic energy to electron thermal energy (e.g., 
Cargill & Papadopoulos 1988), so we neglect these effects. 
Our treatment therefore maximizes the possible offset be- 
tween T c and 2]. 

We work with a flat A-dominated cold dark matter cos- 
mology with matter density flo = 0.3, cosmological con- 
stant Q\ = 0.7 and expansion rate at the present time 
Ho = 70 km s _1 Mpc -1 . We set the baryon density 
Qh = 0.04 and the normalization parameter as = 0.9. 
Our simulation employs 324 3 cold dark matter particles 
and the same number of non-radiative gas particles in a 
cosmological volume 100 ft -1 Mpc on a side. The mass per 
gas particle is 3.26 x 10 s h~ 1 M & , and the nominal gas mass 
resolution is about 10 10 /i~ 1 Af Q for our choice of the num- 
ber of SPH neighbors, N = 32. We note that the CDM 
initial conditions and the large simulation volume allow us 
to correctly model the merging process during the forma- 
tion of filamentary structures and the assembly of massive 
halos which contribute to heating of the IGM as well as 
the intracluster medium (ICM; e.g., Keshet et al. 2003). 

3. RESULTS 

Figure 1 shows the baryon distribution at z = in a 
slice 20 hr 1 Mpc thick (top), the hot component with 
temperature 10 6 < T < 10 7 K (middle), and gas in which 
T c < 0.57]. The latter two components are quite sim- 
ilar, suggesting that much of the hotter WHIM has a 
two-temperature structure. The distribution of the two- 
temperature gas appears quite complex even in clusters, 
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Fig. 2. — Electron temperature plotted against the mean gas 
temperature. The number density of gas particles is shown in grey- 
scale. The rectangular box indicates the so-called WHIM which we 
define to have 10*^ < T c < 10 7 K. 



indicating that shocks induced by mergers during hierar- 
chical assembly propagate and leave two-temperature gas 
within clusters. We have checked the time variation of 
the specific entropy of the gas elements and verified that 
the two-temperature regions have recently been shock- 
heated. In the gas around large-scale structure, continu- 
ous shock-heating disturbs equipartition sufficiently close 
to the shock, while more distant gas is often relaxed. 

Figure 2 shows the electron temperature T c plotted 
against the mean temperature T = (n e T e +mTi)/(n e + ni), 
which is the temperature we would obtain in a single tem- 
perature model. Interestingly, the two-temperature devia- 
tion is largest in the WHIM regime (indicated by the rect- 
angular box). Note that there are some gas elements with 
T e <C T{. This is because our simulation does not include 
electron heating processes other than Coulomb collisions. 

Figure 3 shows the radial distribution of T c (solid line) 
and T (dashed line) for one of the most massive halos in 
our simulation. We used a friends-of-friends group-finder 
to locate dark matter halos and define the virial radius 
such that the mean matter density within this radius is 
200 times the critical density. The halo we consider here 
has a virial mass of 3 x 10 14 ft _1 M Q and a virial radius 
of 1.1ft -1 Mpc. In Figure 3, the difference between the 
two temperatures is large at r > 1ft -1 Mpc and appre- 
ciable even at r < lh^ 1 Mpc. The overall profile looks 
similar to those found in simulations of single clusters by 
Takizawa (1998) and Chieze et al. (1998). Owing to the 
difference in relaxation time (see eq. [3] ) in the dense cen- 
tral part and in the low density outer regions, the electron 
temperature profile appears steeper than that of the mean 
temperature at r > O.l/i -1 Mpc. The bottom panel of 
Figure 3 shows the baryonic mass fraction as a function of 
both T c and T (thick and thin histograms, respectively). 
We selected gas particles within 5 h^ 1 Mpc of the cluster 
center. The distribution of T e is skewed significantly to 
smaller values. The total gas mass contained in this re- 
gion (r < 5k- 1 Mpc) is 1.0 x lO 14 ^ 1 M Q , and nearly 50% 
of it has T c < 1 keV; i.e., the mass of the warm 
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Fig. 3. — Radial distribution of the electron temperature (solid 
line) and the mean temperature (dashed line) in one of the most 
massive clusters in our simulation (top panel). Bottom panel shows 
the mass fraction of gas within 5h~ Mpc from the cluster center in 
temperature bins, using the electron temperature (thick histograms) 
and the mean temperature (thin histograms). 

component surrounding the cluster is greater than the gas 
mass within the cluster itself. 



4. COSMOLOGICAL IMPLICATIONS 

The two-temperature structure of the WHIM and ICM 
has many important implications. For example, the rel- 
ative fractions of metal ions are of considerable interest 
for WHIM observations. Figure 4 shows the fraction of 
Ovi, Ovn, and Ovin ions as a function of temperature 
assuming collisional ionization equilibrium. 4 Clearly the 
abundances have a strong dependence on electron temper- 
ature, because the populations are primarily determined 
by electron collisions (photo-ionization is unimportant for 
these high-level ions unless T <C 10 5 5 K or the incident 
X-ray intensity is high). The bottom panel indicates the 
possible errors caused by assuming a single temperature 
for the WHIM. A factor of two systematic shift in temper- 
ature, typical of the offsets we find between T and T e near 
shocks, can lead to significant over/under-estimates of the 
abundances. 

We further quantify the importance of non-equipartition 
with a model oxygen line emission map. We put the super- 
cluster region indicated in Figure 1 at z = 0.03 and com- 
pute the surface brightness of Ovn emission following the 
procedure of Yoshikawa et al. (2004). We adopt a sim- 
ple relation between the gas metallicity and the local gas 
density of the form Z = 0.02(p gas /p) - 3 Z©, where p is 
the mean baryon density. We then compute the surface 
brightness for both T c and T. Figure 5 shows the 



4 The assumption of collisional ionization equilibrium may not be strictly valid in the WHIM. We nevertheless adopt this commonly-made 
assumption to compute the line cmissivities in a post-processing manner. 
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Fig. 4. — Relative abundance of oxygen ions /q v , (solid), /o v it 
(dashed), and /o V m (dotted) as a function of electron temperature. 
The bottom panel shows the ratio /(T)//(T/2) which indicates the 
error induced by a factor of two shift (misestimate) in T e . 



derived profiles. While the overall results appear sim- 
ilar in the two cases, there are substantial local devia- 
tions. In particular, around some bright spots, the emis- 
sivity is greater by up to an order of magnitude when 
computed using the electron temperature. The ratio of 
the surface brightness S(T C )/S(T) is shown in the bot- 
tom panel of Figure 5. There, we mark the regions with 
S(T C ) > 10~ n ergs -1 cm -2 sr~ 1 , which is the nominal 
detection limit of the proposed DIOS mission (Yoshikawa 
et al. 2003). 

Similar considerations may also apply to efforts to study 
the WHIM using UV observations (e.g., Furlanetto et al. 
2004), although the steep temperature dependence of i e i 
implies that equipartition is more accurate in this regime 
(see Fig. 2). In principle, one can measure the degree of 
equilibration using various line intensities from far UV to 
soft X-rays (e.g., Ghavamian et al. 2001; Raymond et al. 
2003) as well as the metal line Doppler widths. Measuring 
the age of the WHIM from the degree of equilibration will 
put a stringent constraint on theoretical models of large- 
scale baryonic structure formation. 




Fig. 5. — Surface brightness of Ovu emission for the region in- 
dicated in the middle panel of Figure 1 using T e (solid line) and T 
(dot-dashed line). Bottom panel shows the ratio S(T e )/S(T). Dots 
in the bottom panel indicate the region where S(T e ) is greater than 
10~ n ergs -1 cm -2 sr — x . 

Although our simulation does not include radiative cool- 
ing, we note that the lower electron temperature can have 
a significant impact on the radiative cooling efficiency in 
dense regions. The gas cooling rate has a very steep tem- 
perature dependence at 10 5 < T < 10 6 K regardless of 
the metallicity (Sutherland & Dopita 1993). As shown in 
Figure 2, recently shock-heated gas with T ~ 10 6 K may 
have an electron temperature T c ~ 10 5 K. Such gas should 
cool rapidly and in a highly non-equilibrium manner. This 
warrants further study of the effect of the two-temperature 
structure on the IGM using a simulation with radiative 
cooling. 

Overall, our simulation suggests that theoretical studies 
of the WHIM and correct interpretation of the observa- 
tional data from current and future X-ray missions requires 
explicit consideration of these relaxation processes. 

We thank Kohji Yoshikawa for providing emissivity 
data. The simulations were performed at the Center 
for Parallel Astrophysical Computing at the Harvard- 
Smithsonian Center for Astrophysics. 
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